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RARE EARTH OPTICAL TEMPERATURE SENSOR 


Donald L. Chubb and David S. Wolford 
National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

A new optical temperature sensor suitable for high temperatures (>1700 K) and harsh environments is 
introduced. The key component of the sensor is the rare earth material contained at the end of a sensor that is in 
contact with the sample being measured. The measured narrow wavelength band emission from the rare earth is 
used to deduce the sample temperature. A simplified relation between the temperature and measured radiation was 
verified experimentally. The upper temperature limit of the sensor is determined by material limits to be 
approximately 2000 °C. The lower limit, determined by the minimum detectable radiation, is found to be 
approximately 700 K. At high temperatures 1 K resolution is predicted. Also, millisecond response times are 
calculated. 


I. INTRODUCTION 

There are a limited number of temperature sensors suitable for high temperatures (>1500 °C) and harsh 
environments. Platinum-rhodium type thermocouples, 1 which can operate in reactive environments, are suitable for 
temperatures up to 1700 °C. However, similar to all thermocouples, they are not suitable for electrically hostile- 
environments. For temperatures beyond 1700 °C radiation thermometers 1 are generally used. However, radiation 
thermometers require knowledge of the emissive properties of the sample being measured. For a sample with 
constant emittance (gray body) a radiation thermometer can be used without knowing the emittance. 1 However, for a 
sample with an emittance that depends on the wavelength, a radiation thermometer is not suitable. 

In this paper we report on a new optical temperature sensor suitable for high temperatures and harsh 
environments. This sensor does not require knowledge of the emissive properties of the sample being measured. The 
key to the operation of this sensor is the narrow band emission exhibited by rare earth ions (Re) such as ytterbium 
(Yb) and erbium (Er), in various host materials. Depending on the host material this sensor will operate at 
temperatures greater than 1700 °C. 

Most atoms and molecules at solid state densities emit radiation in a continuous spectrum much like a 
blackbody. However, the rare earths, even at solid state densities, emit radiation in narrow bands much like an 
isolated atom. The reason this occurs is the following. For doubly (Re ++ ) and triply charged (Re +++ ) ions of these 
elements in crystals the orbits of the valence 4f electronics, which accounted for visible and near infrared emission 
and absorption, lie inside the 5s and 5p electron orbits. The 5s and 5p electrons "shield" the 4f valence electronics 
from the surrounding ions in the crystal. As a result, the rare earth ions in the solid state emit in narrow bands, much 
like the radiation from an isolated atom. The rare earths of most interest for the optical temperature sensor have 
emission bands in the near infrared (800 < X <3000 nm). 

Development of the rare earth optical temperature sensor has resulted from research on rare earth 
containing selective emitters for thermophotovoltaic (TPV) energy conversion. 2 In that research we have found that 
rare earth doped yttrium aluminum garnet (Re x Yt3. x Al 5 0]2) where Re=Yb,Er,Tm or Ho is an excellent selective 
emitter. It is chemically stable at high temperatures (>1500 C) and produces emittances of e(X) ~ 0.7 in the emission 
bands. 

In the following section the theory of the operation of the sensor is presented. Following that discussion, 
experimental results verifying the sensor operation will be presented. In the final section conclusions will be drawn. 


n. THEORY OF OPERATION 

Figure 1 is a schematic drawing for the rare earth optical temperature sensor. The sensor consists of 
5 components; a rare earth containing end piece, an optical fiber, a narrow band optical filter, an optical detector, 
and electronics to convert the detector output to a temperature. The rare earth containing end piece, which is in 
contact with the sample to be measured, is attached to the optical fiber. Radiation from an emission band of the rare 
earth, which is proportional to the sample temperature, passes through the optical fiber to the bandpass filter. The 
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narrow band filter transmits to the detector only wavelengths within the emission band of the rare earth. Output from 
the detector is then converted to a temperature by an analog electronics package. 

The upper temperature limit for this sensor will be determined by the temperature limits of the materials. If 
ytterbia (Yb 2 0 3 ), which has a melting point of 2227 °C is used for the emitting end of the optical fiber then 
temperatures «2000 °C should be possible if a yttria (Y 2 0 3 , melting point = 2410 °C) optical fiber is used. If a 
sapphire (A1 2 0 3 ) optical fiber is used then the upper limit is reduced since the melting point of sapphire is 2072 °C. 
These materials are chemically stable in atmosphere at high temperature. 


A. Relation Between Temperature and Detector Output 


Assuming the rare earth containing end piece emits uniformly over the area of the fiber, A,, the radiation 
power, q 3 (^,T 8 ) entering the optical fiber from the rare earth at wavelength, X, is the following. 


qs(^ T s) = e b (A)e B (A,T s )A t w/nm 


( 1 ) 


Where E h (X) is the hemispherical spectral emittance into the optical fiber and e B (>.,T s ) is the blackbody 
emissive power 3 at the temperature, T s , of the rare earth containing end piece. 


e B (*.T s ) = 


27TCi 


\ 5 [exp(c 2 t X.T S ) — l] 


W / nm ■ cm ‘ 


( 2 ) 


cj =h C Q =5.9544 xlO 15 Wnm 4 /cm 2 


C 2 = hc 0 /k = 14.388 xl(r nmK 


(3a) 

(3b) 


Where c Q is the speed of light in vacuum, h is Plank's constant and k is Boltzmann’s constant. 

We assume that the emission band of the rare earth is wider than the bandwidth (X^ < X < ^ ( ) of the optical 
filter. In other words, e b (X) > 0 for X^ < X < X,. Also, the transmittance of the optical filter, Xf(?i) = 0 for > X > X t 
Therefore, if the optical filter transmittance is Tf(X), then the power impinging on the detector, q d (A,,T s ), is the 
following. 


^d^’Ts) " T f (^) T £ (^) e b(^) e B(^T s )A/Ff d ^X<X{ (4a) 

q d (X,T s ) = 0 \ U >X>\ ( (4b) 

The term F d is the fraction (F, d < 1) of the radiation that leaves the end of the optical fiber and reaches the detector. 

It is a geometrical factor that does not depend on X or T s and is called the view factor or configuration factor 3 for the 
fiber to the detector. Included in t f (X) is the loss of radiation that escapes out the sides of the optical fiber, as well 
as, absorption losses. 

In using equation (4) for the radiation power arriving at the detector the following approximations are being 

made. 


(1) radiation leaving rare earth containing end piece is uniform across A,. 

(2) The filter transmittance, Tf, optical fiber transmittance, x ( , and the spectral emittance, e b , depend only 
on wavelength. 

(3) Only radiation originating from the rare earth end piece reaches the detector. 

If the rare earth containing end piece is of uniform thickness and at a uniform temperature in the direction parallel to 
the surface of the end piece then approximation 1 will be applicable. Since q fi (A.,T S ) is small absorption in the optical 
filter and fiber wifi be small so that Tf and x ( will be independent of q s (X, T s ) and depend only on X. However, if a 
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significant temperature drop (AT/T S > .05) occurs across the rare earth containing end piece then e b will be a 
function of T s , as well as, X 2 To minimize this effect the rare earth containing end piece must be thin (<0.05 cm), 
but not too thin. The emittance, £*,, depends on the film thickness. 2 If the film is too thin (<0.005 cm) the emittance 
will be reduced to a value such that q s is too small to be detected. The indices of refraction of the optical fiber and 
rare earth containing end piece will change slightly with temperature. We neglect this effect on e b and T ( . 

If the sample being measured has large emittance within the emittance band of the rare earth then this 
sample radiation will contribute to q s . This background radiation can be eliminated by placing a low emittance 
material such as platinum (Pt) between the sample and the rare earth containing end piece. In the experiment to be 
discussed in the next section, Pt foil was placed between the sample and the rare earth containing end piece. Any 
radiation, that impinges on the sides of the optical fiber will not reach the detector. This radiation will merely pass 
through the fiber in wavelength regions where the fiber is transparent or be absorbed in wavelength regions of high 
absorption (X > 5000 nm for sapphire). None of this radiation will be refracted such that it can reach the detector as 
long as the index of refraction of the fiber is greater than the index of refraction of the surroundings. 

Now consider the relation between the detector output and the temperature, T s . A photovoltaic (PV) 
detector, such as silicon, is used to convert the radiation input into a current output. A photoconductive detector 
could also be used as the optical detector. However, a silicon PV detector was chosen because of its superior 
detectivity, time constant and low cost. 4 The equivalent circuit 5 for a photovoltaic device is shown in figure 2. The 
current generated by the input radiation is i ph , the current that flows if a potential V is applied when the device is in 
the dark is i dark and the output current that flows through the load is i. The series resistance is R s and the shunt 
resistance is R sh . Applying Kirchhoff s Law to the circuit loop containing R s and Rg h yields the following. 

R sh(*ph ” i dark ) — 

i = ( 

R sh "h R s 


The dark current is the following. 5 


idark =i sat[ ex P( aV / T d)- 1 ] a - ^ = 8.62x 10 5 (6) 

Where i sat is the so-called saturation current, e is the electric charge and T d is the temperature of the device. Since 
V = V L + iR s equation (5) becomes the following. 

i= - 1 {Rsh [iph - i sat (expfaV L /T d ) exp(ai R s /T d ) - 1 )] - V L } (7) 

Since this current output will be feed into an operational amplifier that acts like a short circuit (V L = 0), the 
current is the following. 


i S c = P ~T~ { i P h - 1 sat [ ex p(ai sc R s /Td ) - 1]} (8) 

K sh + K s 

We now assume that the photovoltaic detector has very small series resistance (R* —» 0). This is a good 
approximation for silicon detectors. 5 In that case the exponential term in equation (8) goes to one. Therefore, the 
short circuit output current equals the photon generated current. 

The photon generated current is given by the following equation. 5 

i sc = i ph — Sr(X)q d (A.,T s )dA. (10) 

Appearing in equation (10) is the spectral response of the PV detector, Sr(X), which gives the current generated per 
unit of radiation power incident on the PV detector. In using equation (10) we are assuming that the PV detector 
responds to all radiation within the limits (A^ < X < X ( ) of the optical filter. In other words, St(X) > 0 for X u < X < X t . 
Substituting equation (4a) in equation (10) yields the following result for i*.. 
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i sc (T s ) = A,F, d f N f (X)t , (X* b (X)Sr(X)e B (X, T s )<& 


( 11 ) 


To proceed further we must have information on the wavelength dependence of if, T,, E b and Sr. The most 
simplification occurs if the bandwidth of the optical filter is so small that T f , T,, E b , Sr and e B are constant for 
A* < A < X , . In that case equation (11) becomes the following. 


1 sc (T s )- A/F, d x f T(EbSr 


2tic j 

X / - x u 

_exp(c 2 /X f T s )-l_ 

* 5 f . 


( 12 ) 


Where Af = (Xu + A)/2 is the center wavelength of the optical filter. For most all applications Af < 2000 nm and 
T s < 2500 K. Therefore, c 2 /AfT s > 2.88 and exp[c 2 /A f T B ]>>I. Also, A^F^ = A ( F, d 3 , where Ad is the detector area and 
F d is the detector to optical fiber view factor. In that case equation (12) becomes the following. 


1 sc (T§ ) — 2 Tic i A d F d / ii f T f E b Sr 



exp[-c 2 A f T s ] 


(13) 


Solving this equation for 1/T S yields the following result. 

± = h. [1„ c -toiler,)] 

4 c 2 


(14) 


Where, 


C = 


2rtc,A d F d( T f T f e b Sr 



(15) 


The constant C, which is independent of T„ can be determined by a calibration procedure. At some known 
calibration temperature, T c , the short circuited current i sc (T c ) is measured. Therefore, 


ln C = r^r + lni sc (T c ) 

A f 4 


(16) 


and 


T = 

1 s 


~ + ~ {l n [l sc (Tc)]- ln [i S c (T s )]} 


-1 


(17) 


Therefore, with the appropriate analog electronics the measured short circuit current, i sc , of the PV detector can be 
converted to a temperature, T s . 

Remember that equation (17) was derived assuming that Tf, x ( , E b and Sr are constants for A^ < A < A,. This 
is a good approximation for the optical fiber transmittance, T, . The reason being that the most probable optical fibers 
for high temperature are sapphire or yttria (Yt 2 0 3 ) which have nearly constant transmittance for wavelengths of 

interest (800 < A < 2000 nm). It is not a good approximation to assume Tf, E b and Sr are independent of X. However, 
as shown in appendix A where wavelength dependence of Tf, E b and Sr is included the short circuit current, i 8C , can 
be closely approximated by the following expression. 
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i sc = C 0 exp[-c 2 /A, f T s ] 


(18) 


Where C 0 is independent of T s and depends on x^A), e b (A) and Sr(A), where Af is the center wavelength of the optical 
filter. As a result, equation (17) applies even when wavelength dependence of T f , Eb and Sr are included. If the 
bandwidth, AA f , of the optical filter is too large (AAf > 25 nm) then C 0 will depend on T s . Therefore, making AAf 
small results in C 0 being nearly independent of T s . For a hypothetical sensor that uses a ytterbium containing emitter, 
which has an emission band centered at A E = 955 nm, and an optical filter with Af = 950 nm and AAf =10 nm the 
parameter C 0 varies less than 0.65 percent for 620 < T s < 2500 K. If the filter bandwidth AAf = 25 nm then C 0 varies 
less than 3.0 percent for the same temperature range. As the results in Appendix A show the variation in C 0 for high 
temperatures (1220 < T s < 2500 K) is much less. For the AAf = 10 nm filter the variation is less than 0.1 percent and 
for the AAf = 25 nm filter the variation is less than 0.5 percent. 


B. Temperature Error 

Consider the temperature error that results from using equation (18) for the short circuit current. Let i sc t(T st ) 
be the short circuit current corresponding to the true temperature, T sl . Let i SC e(T S e) be the short circuit current that has 
a corresponding temperature, T se , and is in error from the true short circuit current, i scl , by the factor, p. As a result 
the following expression applies. 


* see (^se ) — pi set (Tst ) ^ 

Using equation (19) and equation (17) the following result is obtained for the error in temperature, (T st - T se )/T st . 


T — T 

1 st 1 se 


— T st In p 
c 2 

— T st lnp-1 
c 2 


And since V c 2 T st In p «1, 


T st — T se L 

St - " »■ — T si ln P 

T st c 2 


( 20 ) 


In the case of the AA f = 25 nm and Af = 950 nm optical filter just discussed, the maximum error in i^ that results 
from using equation (18) is less than 3 percent for 620 < T s < 2500 K. Therefore, p = 0.97 and 


T st T se < 2.0 x 10~ 6 T c , 


620 <T S <2500 K 
A f = 950 nm, AA f =25 nm 


( 21 ) 


Therefore, if T st = 2500 K then T st - T se < 12.6 K is the maximum possible error in T s for AAf = 25 nm optical filter 
for 620 < T s < 2500 K. For the AA f = 10 nm and Af = 950 nm optical filter the maximum error in i sc is less than 
0.65 percent for 620 < T s < 2500 K. Therefore, p = 0.9935 and 



<4.3xlO -7 T st 


620 <T S <2500 K 
Xf =950nm, AX f = 10 nm 


( 22 ) 


For T s , = 2500 K the maximum possible error in T s is T s - T se < 2.7 K. If we restrict the temperature range to 
1220 < T, < 2500 K then p = 0.995 for the AAf = 25 nm filter and p = 0.999 for the AA, = 10 nm Filter. As a result the 
following results are obtained. 
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T st -T se < 2.1K 


kf = 950 nm, AXf = 25 nm 


T sl = 2500K 


(23a) 


Tst ~T se < 


Xf = 950 nm, AXf = 10 nm 


T st = 2500K 


(23b) 


From these results we conclude that temperature errors resulting from the use of equation ( 1 8) for i sc will be less 
than IK for a temperature sensor that uses the Xf = 950 nm, AXf = 10 nm Filter over most of the useful temperature 
range. 


C. Temperature Limits 

What is the useful temperature range of the sensor? The upper temperature limit is set by the durability of 
the material and has already been discussed. However, the lower temperature limit is set by the minimum short 
circuit current, i^, that can be measured. To estimate this lower temperature limit we make use of the results of 
appendix A. Assume that the optical filter and Yb containing emitter have the following characteristics, Xf = 950 nm, 
AXf = 10 nm, XfMAx = 0.35, X E = 955 nm, e MAX = 0.8, s' = 0.01. Also assume a silicon detector of 1 mm in diameter 
(r d = 0.05 cm) and maximum spectral response Sr(X g ) = 0.5 AAV, where Xg = 1 100 nm is the wavelength that 
corresponds to the silicon bandgap energy, E g = l.leV. In addition, assume the optical fiber transmittance is x, = 0.2 
and the optical fiber diameter is approximately the same as the detector diameter and is as close as possible to the 
detector so that the detector to optical fiber view factor is F d ~ 0.5. In that case the results from appendix A yield 
C 0 - 5.8 A. The minimum value i sc that can be detected will be determined by the electronics package and the 
electronic noise in the system. Using operational amplifiers 6 in the electronics package should allow minimum short 
circuit currents of i sc = 10 8 A to be easily detected. In that case with C 0 = 5.8 A and Xf = 950 nm equation (18) yields 
T s = 750 K. Obviously, at a given temperature a larger value of AXf will result in larger C 0 and therefore larger i K . 
Thus, for AXf = 25 nm with all other conditions the same yields C Q = 13.8 A and a temperature limit of T* = 720 K. 
However, it must be remembered that AXf = 25 nm results in a larger error in T s than the AXf = 10 nm case. 

Based on the results just discussed, a lower temperature limit of the order of 700 K should be possible with 
a silicon detector and optical fiber of 1 mm diameter. Using a larger diameter optical fiber and detector will result in 
a lower temperature limit. However, the larger optical fiber will also result in a larger temperature error since 
thermal conduction away from the sample will be larger. 


D. Response Time 

The time behavior of the output of the rare earth containing end piece of the sensor, q s , depends on how fast 
the end piece reaches a steady state temperature after a change in the sample temperature, T s (t). This response time, 
T, will be much larger than the transmission time of q s through the optical fiber and Filter since they occur at the 
speed of light. The response time of the detector and electronics should also be much smaller than x. Therefore, the 
response time of the end piece, x, will determine the response time of the sensor. 

Assume the end piece behaves in a one-dimensional manner so that its temperature is governed by the one- 
dimensional heat conduction equation. 


fcth-T+PCp^—O (24) 

dx z ot 

Where t is time, x is the coordinate perpendicular to cross sectional area, A f , of the end piece, k^ is the thermal 
conductivity, p is the density and c p is the specific heat. If equation (24) is written in dimensionless form using 
dimensionless length x = x/d , where d = thickness of the end piece then we find that the appropriate dimensionless 

time is t = t / x 0 where, 
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( 25 ) 


pc d 2 
x o =—r — 

k th 

In reference 7 solutions to equation (24) are presented for various boundary conditions and initial conditions. In all 
cases the time dependence is the following, where a is a constant. 

T - e “ at/T ° (26) 

This same result can be obtained by using the method of separation of variables on equation (24). The constant a 
depends on the boundary conditions and is the order of 1 . Therefore, equation (25) can be used to estimate the 
response time, x, of the sensor. 

As equation (25) indicates the response time is a quadratic function of d so that by making the end piece 
thin the response time will be small. Assume that the end piece is a rare earth oxide such as ytterbia (Yb 2 0 3 ) or erbia 
(Er 2 0 3 ). Then p < 10 gm/cm 3 and c p < 1 J/gmK. 8 There is no data available on the thermal conductivity, k< h , or the 
rare earth oxides at high temperature. However, we expect they will have low thermal conductivity at high 
temperature similar to other ceramic materials such as alumina where ~ 0.1 W/cmK. Therefore, assume 
^ ^ 0.05 W/cmK so that equation (25) yields the following for the response time, X. 

T<200 d 2 sec din cm (27) 

Thus, if the minimum thickness end piece is d = 0.005 cm then the response time is x < 5 msec. As discussed earlier, 
the response time of the sensor should be determined by the end piece response time, x. 


III. EXPERIMENTAL VERIFICATION OF SENSOR OPERATION 

To verify that equation (17) applies for determining the temperature, T s , the experiment shown in figure 3 
was used. The sensor consists of a 0.055 cm diameter sapphire fiber with an erbium aluminum garnet (Er 3 Al 5 02) 
emitter of thickness, d = 0.03 cm attached with platinum (Pt) foil. Erbium has an emission band at X = 1000 nm. The 
Pt foil serves two purposes. First of all it holds the emitter to the end of the fiber and secondly it blocks all radiation 
from the sample being measured. The thin (=0.005 cm) Pt foil will produce negligible temperature change between 
the sample and the emitter. In this case the sample being measured is Pt foil with a type R platinum versus platinum 
(13 percent) rhodium thermocouple attached behind the Pt foil. The Pt sample is held with an alumina holder, which 
is located in the center of a high temperature atmospheric furnace. At the other end of the sapphire fiber is a chopper 
and monochromator, which serves as the optical filter with a bandwidth of AXf * 2 nm. A silicon detector converts 
the radiation leaving the monochromator to a short circuit current, i^, which is measured with a lock-in amplifier. 

The experimental procedure to verify equation (17) was as follows. First i sc (T c ) was determined by 
measuring i^ when the thermocouple was at temperature T c . Knowing i sc (T c ), T c and Xf = 1012 nm, equation (17) 
was then used to calculate T s for a series of different thermocouple temperatures, T tc . If equation (17) is valid then 
T s = T xc . In figure 4 the sample temperature, T s , calculated from equation (17) is plotted as a function of the 
thermocouple temperature, T tc , for two different ranges. In figure 4(a) the range is 863 < T tc < 1430 K and in 
figure 4(b) the range is 1334 < T tc < 1879 K. 

As figure 4 shows the agreement between T s and T tc is excellent. The largest error, err = |T T c - T s |/T tc , 

occurs at the lowest temperatures. For figure 4(a), err MAX = 0.03 and for figure 4(b), err MAX = 0.008. As a result of 
this close agreement between T s and T tc we conclude that equation (17) is valid for determining the sample 
temperature, T s , from the measured short circuit current, i sc . 

Figure 5 shows the lock-in amplifier output, which is directly proportional to the short circuit current, i sc , as 
a function of T s for the same conditions as figure 4(b). Two things should be noted. First of all, because of the 
exponential dependence of i sc on T s (eq. (18)) the sensor sensitivity increases with temperature. As figure 5 shows 
for high temperatures (>1 100K) a small change in T s results in a several millivolt change in amplifier output. Thus, 
we expect IK temperature resolution at high temperature. The second thing to note is even though the 
monochromator bandwidth is small (=2 nm), the amplifier output is many millivolts. Therefore, the electronic circuit 
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that converts to T s (eq. (17)) will operate with a large input signal, which should eliminate electronic noise 
problems. 


TV. CONCLUSION 

The operation of the rare earth optical temperature sensor has been experimentally verified. For the 
temperature range 863 < T < 1430K the maximum temperature deviation of the sensor temperature from a measured 
thermocouple temperature was 3 percent. For the temperature range 1334 < T < 1879K this error was reduced to 
0.8 percent. Calculations show that a simple relation (eq. (18)) between temperature and detector output can be used 
to determine the temperature. The error resulting from using this relation will be less than IK. 

Material durability limits the upper temperature of the sensor. Using a yttria optical fiber and rare earth 
oxide emitters the upper temperature limit should be approximately 2000 °C. The lower temperature limit is 
determined by the minimum signal that can be detected. This limit is calculated to be approximately 700 K if a 
silicon detector is used. 

Response time of the detector will be determined by how fast the emitter responds to a temperature change. 
This response time was estimated using the equation for thermal conduction. Results indicate that response times the 
order of msec are possible. Because of the exponential dependence of sensor output on the temperature we expect 
temperature resolution of 1 K to be possible. 

The electronics package to convert the sensor output to a temperature is under development. This 
electronics package can be designed using commercially available parts. Thus the most expensive component of the 
sensor is expected to be the optical fiber. 
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APPENDIX A. RELATION BETWEEN TEMPERATURE AND DETECTOR SHORT CIRCUIT CURRENT 


The complete relationship between the sample temperature, T s , and the PV detector short circuit current is 
given by equation (11). 


dX 

i sc (T s ) = 27Cc 1 A,F (d J x f (X)r,(X)e b (X)S r (X) , „ 7 dO 

K X 5 [exp(c 2 /XT s )-l] 

As stated in the discussion about equation (17), it is a good approximation to assume the optical fiber transmittance, 
T,(X), is constant for Xy< X < X,. The spectral response of a PV detector can be approximated as a linear function of 
wavelength. 9 


Sr(X) = -~-Sr K — 
X g S W 


0 < X < X g 


Sr(X) — 0 X g < X 

Where Sr g is the spectral response at X = Xg = hcJEg and E g is the bandgap energy of the PV detector. A linear 
approximation can also be used for the spectral emittance, E b (X), of the rare earth selective emitter 2 for a narrow 
wavelength region around an emission band centered at X = X E that is located within the transmission band 
X [X <X<X i of the optical filter. 


E b(k) - e MAX -X) X u < X < X E 

Cb(^) = e MAX — — Xg ) X E < X <Xf 


Where e max = £(X E ) is the maximum emittance and £' is the slope of the linear £b(X) versus X approximation. For a 
narrow bandwidth interference filter the transmittance, Tf(X), can be closely approximated by a Gaussian function. 


T f (^)- T fMAX ex P' 


-ln2 


2(X-X { ) 


-|2 


AX f 


(A3) 


Where Xf is the center wavelength of the transmission band, TfMAX = Tf(Xf), and AXf is the bandwidth defined by the 
wavelengths X& = Xf - AXf/2 and Xf, = Xf + AXf/2 where Tf (Xf ± AXf/2) = Tf max/2. Figure A1 compares the 
experimental transmittance of an interference filter with Xf = 949 nm and AXf = 14 nm with equation (A3). As can be 
seen the agreement is quite good. 

Figure A2 shows the approximations for Sr(X)/Sr g , Eb(X)/£ M Ax» and Tr(X)/TfMAx given by equations (Al), 
(A2) and (A3) for the case where Xf = 950 nm, AXf = 10 nm, Xe = 955 nm, £ MA x = 0.8, £jJ = 0.01 nm -1 (ytterbium) 
and X g = 1 100 nm (silicon). These approximations were used in equation (11) and i sc was calculated by numerical 
integration. The integration limits were Xu = Xf - 4AXf and X, , = Xf + 4AXf, which insures that at these limits the 
integrand in equation (11) vanishes. To determined if equation (18) is a valid approximation for i*. we calculated C D 
in equation (18). 


c 0 =i sc exp[c 2 /X f T s ] (A4) 

Where i* in equation (A4) is obtained by numerical integration of equation (1 1). 

Figure A3(a) shows C 0 as a function of T s for the AXf =10 nm filter (Xf = 950 nm, TfMAX = 0.35), fiber 
transmittance, X = 0.2 and the ytterbium emitter (X E = 955 nm, e M ax = 0.8, = 0.01 nm' 1 ) and silicon detector 
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(X g =1100 nm, Sr g = 0.5 AAV). As can be seen, C 0 changes only by a small amount for 620 < T s < 2500 K. Thus, C 0 
is nearly independent of T s and equation (18) is a valid approximation for i sc . Figure A3(b) shows C 0 for the case 
where AXf - 25 nm with all other conditions being the same as in figure A3(a). In this case C 0 changes more than the 
AXf = 10 nm case, but equation (18) should still be a good approximation for i sc . If the optical filter and emission 
band are matched (Xf = X E ) then C D varies even less over the temperature range (620 < T s < 2500 K). For the 
AXf =10 nm filter the maximum variation in C Q is less than 0.65 percent if X E = 955 nm and Xf = 950 nm but is only 
0.33 percent if Xf = X E = 955 nm. 
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Figure 1 . — Schematic of rare earth optical temperature sensor. 
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Figure 3. — Experiment to verify rare earth optical temperature sensor operation. 
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Sample temperature, T s> K £ Sample temperature, T, 


* 



Thermocouple temperature, Tjq, K 



(b) Thermocouple temperature, Tjq, K 

Figure 4. — Comparison of sample temperature, T s , 
determined using equation (1 7) for erbium aluminum 
garnet emitter, and sample temperature measured by 
type R thermocouple, Jjq. Center wavelength, 

X f = 1012 nm. (a) 863 < T rc ss 1430 K, T c = 1430 K. 

(b) 1334 < T TC s 1879 K, T c = 1879 K. 
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Sample temperature determined by sensor, 
T s , K 


Figure 5. — Lock-in amplifier output, which is directly 
proportional to short circuit current, i sc , of silicon 
detector, as a function of the measured sample 
temperature, T s , for same conditions as figure 4(a). 
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920 930 940 950 960 970 980 

Wavelength X, nm 

Figure A1 . — Comparison of actual transmittance, Tf, with 
Gaussian approximation of Tf for narrow bandwidth 
interference filter with center wavelength, Xf = 949 nm 
and bandwidth, AXf =14 nm. Actual Tf is solid line and 
Gaussian approximation is dashed line. 
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Wavelength X, nm 


Figure A 2. — Optica! filter transmittance, Tf/Tf max , Yb 
emitter emittance, et/ e max» and silicon detector 
spectral response, Sr/Sr g used to calculate the 
parameter C 0 in eq. (18). 
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5.85 




(b) Temperature, T s , K 

Figure A3. — Parameter in Eq. (18) as a function of 
temperature for the following conditions, r a = .05 cm, 

Fdl = -5, T| = .2, Tf max ss .35, e ma x = -8» = *01 nm - ^, 

Xg = 1 1 00 nm, Sr g = .5 A/W. (a) Filter center wave- 
length, Xf = 950 nm, bandwith, AXf =10 nm. (b) Filter 
center wavelength, Xf = 950 nm, bandwith, AXf = 25 nm. 
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